Long-range attraction between particles in dusty plasma 
and partial surface tension of dusty phase boundary 
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Effective potential of a charged dusty particle moving in homogeneous plasma has a negative part 
that provides attraction between similarly charged dusty particles. A depth of this potential well 
is great enough to ensure both stability of crystal structure of dusty plasma and sizable value of 
surface tension of a boundary surface of dusty region. The latter depends on the orientation of the 
surface relative to the counter-ion flow, namely, it is maximal and positive for the surface normal 
to the flow and minimal and negative for the surface along the flow. For the most cases of dusty 
plasma in a gas discharge, a value of the first of them is more than sufficient to ensure stability of 
lenticular dusty phase void oriented across the counter-ion flow. 
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I. INTRODUCTION 

A gas-solid phase transition observed in different dusty 
laboratory plasmas [1-8]. It counts in favor of presence 
of strong long-range attraction between similarly charged 
dusty particles. By now some physical mechanisms are 
proposed to explain formation of a regular arrangement 
of micron-sized particles embedded in a gas discharged 
plasma. A part of them is based on the account for elec- 
trostatic fields of strata and walls of discharge tubes [9]. 
In doing, so a generality of solution of the problem of 
the effective attractive potential is lost and, in particu- 
lar, the situation of dusty plasma crystal in a thermal 
dusty plasma of gas burner [5] where there are no walls 
and external fields drops out of such interpretations. 

The most promising approach holds that the ion 
streaming motion causes an attractive wake potential 
behind the dust particles. Originally it was developed 
[10, 11] for particular case of supersonic flows which is re- 
alized in the sheath of radiofrequency discharges. Later, 
it was extended with participation of one of the authors 
of Ref. [10] to the case of subsonic ion flows [12], but 
physics of shielding of dust charged particles was sup- 
posed to be strongly modified with mandatory regard to 
anisotropy and asymmetry of the ion temperature in the 
sheath. 

Below is shown that the attraction between simi- 
larly charged particles can be resulted from a dynami- 
cal screening of the Coulomb potential remaining in the 
frame of a single physical mechanism for both supersonic 
and subsonic regimes. The distinction of the dynami- 
cal screening from the static Debye screening is due to a 
motion of particles relative to screening charges (counter- 
ions or electrons). From the physical point of view, this 
effect can be interpreted as a consequence of loss of spher- 
ical symmetry of the Debye screening cloud around the 
moving charged particle resulted in a space charge with 
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the opposite sign forming in its wake. Not only do this 
space charge compensates the particle's potential but 
may also give rise to a local wake potential of the op- 
posite sign. In the supersonic regime the wake potential 
oscillates inasmach of Cherenkov wave generation. 

This effect is well known in the usual electron-ion 
plasma [13-15]. When the charge is moving, the static 
Debye screening modifies so that the potential in the 
oncoming flow of electrons grows approaching to the 
Coulomb potential with an increasing of velocity while 
the potential in the outward flow decreases up to alter- 
nating of its sign. 

Similar effect was found [15, 16] in a system of grav- 
itating masses where the static screening is absent and 
accounting for the dynamical screening gives rise to an 
alternating potential of gravitational interaction. 

The aim of the present paper is to propose the most 
general model of dynamical screening of field of dusty 
particle charge to explain the observable interparticle at- 
traction. It may be supposed that the distance from the 
particle to the attractive minimum of the wake poten- 
tial determines a period of the resulted lattice of dusty 
particles. 

Besides, the interparticle attraction is to give rise to 
a surface tension of a dusty phase interface. Really, this 
explains a sharp non-diffusion character of a particle den- 
sity variation observed at the interfaces in the experi- 
ments on dusty plasma crystals. The most characteristic 
phenomenon of such kind are voids in a homogeneous 
dusty plasma. The characteristic lenticular form of the 
voids may be explained by strong peculiar dependence of 
the surface tension on the surface orientation relative to 
the ion flow. 

The notion of the surface tension of a dusty phase was 
mentioned briefly by Tsitovich [17] but his estimation of 
its value was erroneous. Conceivably that might be the 
reason why the surface tension was not mentioned in the 
subsequent papers by Tsitovich (as a co-author) on the 
theory [18, 19] of spherical voids. 

The outline of the paper is the following. In Sec. II, the 
model for effective non-potential attractive forces of in- 
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teractions between dust particles will be developed. This 
model can explain an appearance of the lattice structure 
both in gas discharge dusty plasma and in thermal dusty 
plasma of gas burner. In Sec. Ill, a surface tension of 
a boundary of dusty phase is estimated and conditions 
of stability of a lenticular void in the dusty plasma are 
examined. 



II. EFFECTIVE INTERPARTICLE 
INTERACTION 

The system under consideration consists of negatively 
charged dusty particles of the charge Q = —Z p e and 
concentration n Pl and positively charged counter-ions 
of the charge e (for definiteness sake we will consider 
singly charged ions) and mass m. For simplicity, neutral 
molecules of a buffer gas are neglected here and dusty 
particles are considered as point charges. 

While forming of the negative charges of dusty parti- 
cles the electrons are condensed on particles' surfaces and 
their concentration in plasma decreases sufficiently. As 
a result, in the case of great Z p n p their density becomes 
so little that the electron Debye radius of screening of 
the ion charge appears to be much greater the ion Debye 
radius of screening. Thus, we can restrict our considera- 
tion to the simplified model of one-component ion plasma 
in which the negatively charged dusty particles are im- 
mersed. 

If the negative test point charge Q is moving in the 
system of ions and dusty particles with the velocity u it 
gives rise to some perturbations of the system state. Due 
to great difference in masses and concentrations of ions 
and dusty particles we can restrict our consideration to 
a perturbation of the ion component only (just as a per- 
turbation of an electron component is taken into account 
only, as a rule, in the problem of screening of an ion 
charge in the electron- ion plasma). Such perturbation 
of ionic subsystem is described by the set of the Vlasov 
equation for a distribution function over coordinates and 
velocities of ions and Poisson equation for an effective 
potential induced by the perturbation of ion density and 
moving test charge. Interactions between ions and neu- 
tral molecules of the buffer gas are neglected here [20]. 

Such simplified model is well studied in the test- 
particle approach to the electron-ion plasma theory [13- 



15] so we can omit intermediate calculations and write 
down at once the resulted form of the effective potential 
of the moving test charge 



. 3l; exp{k- (r-ut)} 
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is the dynamical permittivity of the ion subsystem. Here 
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FIG. 1: The reduced potential (p — $/(Qk) of the negatively 
charged dusty particle at Z = moving relative ions along 
axis Z with the reduced velocity M = u/v = 3. 



v = (2/csTi/m) 1 / 2 is the mean heat velocity of the ions 
and 



W(t) = 1 - erfi(i) + iy/^te' 



(3) 



Let us choose an axis z along u and introduce dimen- 
sionless variables Z = (z — ut)n, X = xk, K = k/n and 
the Mach number M — u/v relative to the ion heat ve- 
locity. Then, in the test particle's frame accounting for 
the cylindrical symmetry, we get 



J 



$(r) = y d(f>J^deJ dKK 2 sine 



exp{iKA} 
K 2 + W(Mcos9) 



Q K ip{X,Z;M), 



(4) 



where 



y{X,Z-M) 



(X 2 + Z 2 )^ 



--I(X,Z;M), 



I(X, Z; M) = ^ J d<j> d6 J dKW{M cos 0) sin 



exp{iKA} 
K 2 + W{Mcos0Y 



(5) 
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A = X sin 6 cos <b + Z cos 0. 



This integral determines the departure of the effective potential from the Coulomb potential to which the first member 
of ip(X, Z; M) corresponds. When the Cauchy integral over K is taken, we find 



I(X,Z;M)=^- 2 



de sin evw 



2i smh(Ay/-W)Ci{-i AVw) + 



+cosh( AVW) (ir + 2i Shi(AVWf) 



(6) 



It is not difficult to calculate numerically the value of 
this integral in an arbitrary range of the parameter M 
variations at any point X, Z. 

Corresponding results of numerical calculations for 
M = 3 are illustrated in Fig.l. It is seen here that there is 
well defined negative minimum of the potential along the 
axis Z at a distance of order of five Debye radius, which 
has been formed as a result of deformation of screening 
Debye cloud in the vicinity of the moving particle. 

In gas discharges, where forming of crystal structure 
of dusty particles was observed, a diffusion ion velocity 
u is determined by their mobility and an electric field 
strength. To define the mobility, we are to take into 
account a tendency of ions to unite with molecules and 
atoms into complexes of the types N± , Of, Ne^ ', Hef,, 
and recharge effect also. According to [22], the diffu- 
sion velocity of the complex Ne^ is of order of 50 mc -1 
in a characteristic for glow-discharge field E/p — 
1 Vcm _1 torr _1 and the heat velocity v = 400mc~ 1 (at 
T = 300 K). Then M = 0.125. For the ion Ne+, the dif- 
fusion velocity decreases near-threefold as a consequence 
of the resonant recharge, then M ~ 0.05. In the case of 
strong field, the ion drift velocity may not only approach 
to the heat velocity but surpass it [22]. Because of this, 
a wide range of values of the parameter M — u/v will be 
considered below (asymptotic case M»l was discussed 
in [10, 11]). 

Observed quasi-crystal structure in a thermal plasma, 
that is, in the dusty flame of gas burner [5, 8, 23], can 
be treated also as a consequence of movement of charged 
particles relative to ions of the flame. To obtain quan- 
titative estimations of the effect, the additional data on 
a gap of velocities of particle and ion flows are neces- 
sary. This gap is determined by conditions of injection 
of the dusty particles into the flame, and in the steady 
state limit it is likely to go to a sedimentation velocity of 
charged particles in the rising flow. 

The graphs of tp(X = 0, Z < 0; M) for different M 
in the wake of charged particle are illustrated in Figs. 
2,3. As can be seen, the depth of the negative minimum 
|<^ TO m| increases with an increasing of M from 0.01 to I 
(see Fig. 4), and its location Z m i n shifts to the particle. 
At low M \Z m i n \ depends on M as \Z m i n \ = 1.6— 1.71nM 
(see Fig. 5). When further increasing of M, the potential 
becomes oscillating and its first minimum moves away 
from the particle. For M > 2 the graph in Fig. 5 
can be well approximated by the expression \Z m i n \ = 



1.8VM 2 — 1, which is close to the experimental data [1- 
4] on the period L of a crystal lattice in a dusty plasma 
if it is identified as L ~ \Z min \. It should be noted here 
that experimental detection of the minimum of L(M) 
corresponding to the minimum of the graph in Fig. 5 
would be a good evidence in favor of the approach under 
discussion. 

The values M > 1 correspond to a supersonic move- 
ment when Cherenkov radiation of ion-acoustic waves 
takes a place that explains an appearance of space os- 
cillations of the wake potential. It should be noted that 
in the range of values of M from 1 to <~ 3.3 an increasing 
of the depth of the first minimum takes a place but at fur- 
ther gain of M the depth of the first minimum decreases 
(see Fig. 4). 

In the incoming flow (Z > 0) the effective potential 
varies continuously with a gain of M from the Debye 
potential ~ exp{— Z}/z to Coulomb one ~ 1/z. This 
indicates that the forward part of the Debye screening 
cloud is not pressed by the incoming flow to the particle 
but blew away from it. 

The negative minimum of the wake potential points to 
presence of a space charge of the opposite sign. Thus, 
the large charged particles moving relative ion subsys- 
tem together attendant space charges constitute a sys- 
tem of dipoles oriented along the direction of the relative 
motion. This gives rise to formation of hierarchies of 
particles-dipoles along lines of the ion flow. 

Moreover, this interaction between dust particles is 
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FIG. 2: The same as in Fig. 1 along the axis Z at M 
0.1, 0.2, 0.3, correspondingly. 
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asymmetric in such a way that attractive force is commu- 
nicated only downstream the ion flow. This situation was 
clearly demonstrated in the experiments [24, 25] where 
upper or lower particle was pushed by the laser beam 
and it was just lower dust particle that fitted its position 
when the upper particle was shifted and not vice versa. 
So heavily non-conservative character of the interparti- 
cle interactions excludes any possibility to introduce an 
effective potential of the interaction. 

Besides, as it is seen in Fig.f , the resulted negative po- 
tential is long-range in the transverse direction X, also, 
that provides an attraction between particles from neigh- 
bor hierarchies and, in its turn, causes the particles from 
neighbor hierarchies to shift relative one another along 
the ion flow for a distance of order of a half dipole length. 
Because of this interaction is much weaker then along the 
flow, the situations can take place when the lengthwise 
attraction appears quite strong for the hierarchies forma- 
tion but too weak for transverse ordering [8] . 

Thus, the movement of the charged particle relative 
ion flow gives rise to anisotropic potential field in which 
an energy of another similar particle of the charge Q is 
U(r) = Q<f>(r) = Q 2 Kip(r). Then, the effective force of 
interaction between the particles is 



F(r) 



dUjr) 
dr 



-Q 2 k 



dr 



(7) 



For the formation of the crystal structure, it is nec- 
essary that the maximum deep of the attraction energy 
\U m in\ = \Q®min\ = Q 2 n\ip m in\ surpass the heat motion 
energy, that is, 



(8) 



where T p is the kinetic temperature of the particle. Ac- 
cording both experimental data [26] and theoretical es- 
timations [27] T p may sufficiently surpass the ion tem- 
perature. For definiteness sake, we take T p = 1 000 , 
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FIG. 3: The same as in Fig. 1 along the axis Z at M 
1, 2, 3, correspondingly. 
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FIG. 4: Dependence of the depth of negative minimum of the 
reduced potential on M. 



30 



§ 25 
1 20 



S 15 
o 

8 10 
a 



10 12 



14 



M 



FIG. 5: Dependence of the position of negative minimum of 
the reduced potential on M. 
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400 cm , then the condition of the hier- 



archy stability becomes |<£ m i n | > 10 . According to the 
approach under discussion (see the graph in Fig. 4), this 
condition is fulfilled when M > 0.01. 

The long-range attraction between likely charged dusty 
particles resting in an equilibrium plasma (M = 0) was 
found in [28]. It is possible that this kind of attraction 
can give rise to forming of crystal structure but it is ir- 
relevant to the observed [1-4] dependence of the lattice 
period on the ion flow velocity and forming (at small M) 
of hierarchies along the flow non-interacting with which 
other. 



III. SURFACE TENSION OF DUSTY PHASE 

One of characteristic feature of a dispersed phase in 
dusty plasma is presence of sharp (non-diffusion) bound- 
ary surface of a dusty cloud ("dusty drop"). Consider- 
able recent attention has been focused on the problem 
of formation of voids ("dusty bubbles") of the dispersed 
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phase. It is reasonable to suppose that both phenomena 
are of the same nature related with presence of a surface 
tension resulted from the attractive interaction between 
charged particles. 

In general, the surface tension coefficient of an inter- 
phase boundary is defined as the integral of a difference 
between transverse p± and longwise p|| pressures over an 
interphase layer of thickness I 

rl/2 

7 = / 0-lOi) -P\\(zi))dz 1 (9) 

J -1/2 

where the axis z\ directed along the normal to the layer. 

On the molecular level, a pressure can be written in 
the form of a virial equation of state and expressed via 
a number density n(zi), virial of a force of intermolecu- 
lar interaction F(r) and radial correlation function g(r), 
where r = r 2 ri is the distance between particles. For 
a bulk equilibrium fluid the function g(r) is calculated 
as a rule with the use of the Born-Green equation. At 
an interphase surface the situation is much more com- 
plex; this being so, notably rough approximations are 
used here. One of them is Fowler's step approximation, 
when the number density n(zi) is regarded as constant no 
within the liquid phase and zero outside the liquid, and 
radial function g(r) within the surface layer is assumed 
to be the same as in the bulk liquid. Then the surface 
tension coefficient is determined by the Fowler formula 
[29, 30] 

7 = -— o / r^F^g^dr, z > 0. (10) 
° Jo 

It is implied here that the z-component of the vector r is 
normal to the surface. Kirkwood and Buff [31] estimated 
the surface tension of liquid argon on the base of this 
formula and obtained quite satisfactory value (error ~ 
25%). 

It is supposed in the Fowler' formula that the inter- 
molecular interaction is spherically symmetric that does 
not allow to apply this formula as such to the boundary 
of dusty phase. 

Rederiving the Fowler formula with due regard for the 
anisotropy of interparticle interaction we obtain 

7 = -y/ [z 2 F z (v)-xzF x (r)]g(v)dv, z > 0. (11) 

where F z and F x are z and x components of the force 
(7). 

When applying to the boundary of dusty phase, the 
pressures pj_ and py arc to be considered as partial pres- 
sures of dusty component. Then, the surface tension of 
dusty phase should be interpreted as partial property too. 

To estimate its value some simplification will be intro- 
duced. Firstly, the radial distribution will be used in the 
form of the Boltzmann factor 

ff (r)=exp{- — -}, 




FIG. 6: The typical void in dusty plasma. (Copied from 
the Web page http://www.mpe.mpg.de/www_th/plasma- 
crystal/PKE/PKE-Nefedov.html dedicated to the memory of 
A. P. Nefedov. Authorized by Max Planck Institute of Ex- 
traterrestrial Physics (Germany)) 

that is typical for dilute systems. Secondly, when esti- 
mating the surface tension, the interactions within hi- 
erarchies will be regarded only as it is much stronger 
interactions between particles from different hierarchies. 
Then, for surfaces oriented across and along the ion flow 
we get 

7± = - y / z 2 F z (r)g(r)dr, z > 0. (12) 



7|| = Y 7 xzF x (v)g(r)dr, z > 0. (13) 

The axis z in both equations directed along the normal 
to the surface and the forces F z (r) in the equation for 
7_i_ and F x (r) in the equation for 711 are determined by 
interactions along hierarchies. 

A numerical estimation of these expressions for n p ~ 
10 6 cuT 3 and M = 0.01 leads to 71. ~ lOMincm" 1 
and 7|| <~ — 10 4 dincm _1 . Both of these coefficients 
in absolute magnitude surpass sufficiently the typical 
surface tension coefficient of a liquid at normal con- 
ditions jii q ~ 10 2 dincm _1 (as an example, for water 
luq = 70dincm _1 ). Such great difference is quite nat- 
ural. Really, the negative minimum U m of the potential 
U at given conditions is of order of — 10~ 12 erg while the 
depth of the potential well of, for example, the Lennard- 
Jones potential for noble gases is of order of — 10~ 14 erg. 
Inasmuch as the potential U enters (with sign minus) 
into the exponent of the radial distribution function, the 
surface tension coefficient appears to be very sensitive to 
its value. As a result, we have |7j_ i || | 3> ju q in spite of 
7 cx n 2 and n p <C n liq . 

Hydrodynamic stability of a spherical gas bubble or 
liquid drop is determined by value of a Weber number 
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We — pv 2 a/~f, which is the ratio of a dynamic head 
pv 2 to a surface tension pressure 7 /a, where a is the 
radius of the bubble or drop. For the case of the void 
in dusty plasma, the dynamic head pv 2 is determined 
by the flow of ions and neutral gases with the velocity 
v = u = vM~5 - 10 4 M cm c _1 . We can put 7 ~ 7^ for 
the slightly curved surface of a lenticular void with the 
radius of curvature a ~ 2 cm sufficiently greater its size. 
Then at characteristic pressure of order of 0.5torr we have 
We ~ 10 4 M 2 /7j_. The coefficient 7^ increases with in- 
creasing of M, so that the Weber number We attains 
its the most great value at the minimal from considered 
here values of M, i.e. at M = 0.01 when j± ~ 10 5 . 
Then We m ax — 10~ 4 . When the Weber number is so 
small the surfaces of the lenticular void are stable, cer- 
tainly. Their deviations from a spherical form can be due 
to variations of 7 relative to 7j_ as a consequence of in- 
evitable (although small) distortions of the condition of 
orthogonality to the flow on the curved surface. 

Stability of the lenticular void ensures possibility of its 
existence. The problem of its appearance remains to be 
solved. It may be suggested that a phase transition of the 
first kind takes a place in a homogeneous dusty plasma, 
and strong anisotropic surface tension can give rise to a 
peculiar nucleation process. 

IV. CONCLUSIONS AND DISCUSSIONS 

Attractive interactions between likely charged colloidal 
or dust particles in plasma have been discussed in this ef- 
fort in a view to explain the observed crystal structure 
formed from these particles. The crucial point of the 
model under discussion is the presence of the counter- 
ion flow. Its velocity u relative the dust particles is not 
necessary to exceed the heat velocity v of counter-ions 
(or ion-acoustic velocity) as it was supposed in ref. [10]. 
Even at small u when u/v ~ 10~ 2 , the effective attrac- 
tive interaction appears to be sufficiently strong to ensure 
stability of regular hierarchies of particles oriented along 



the ion flow. 

Another consequence of strong attractive interparticle 
interactions is the partial surface tension 7 of the dusty 
phase boundary. 

This concept was briefly discussed by Tsytovich [17] 
who estimated 7 as a work necessary for construction 
of a bulk liquid of the unit surface and height h. As 
a result, he obtained 7 = U m n p h, then he took |J7 TO | = 
100 eV that correspond to our result at M = 1 and found 
7 ~ 10 _2 dincm _1 neglecting the negative sign of U m . If 
he accounted for U m < he would get 7 < 0, that is 
natural as the work for construction of the coupled state 
is to be negative. 

Here, the surface tension is estimated on the base of 
a generalization of the Fowler's formulae taking into ac- 
count the anisotropy of interparticle interactions. As a 
result, the strong dependence of 7 on orientation of the 
dusty interface is found. In particular, when u/v ~ 10~ 2 
we get 7j_ — 10 5 dincm _1 and 711 ~ — 10 4 dincm _1 for 
surfaces oriented across and along the ion flow, corre- 
spondingly. So great positive and negative value results 
in the characteristic lenticular form of the voids in dusty 
plasma having no surfaces along the ion flow. Such a 
form was observed both under micro-gravity condition 
(see Fig. 6) and in terrestrial experiments by Samsonov 
and Goree (see Fig. 5d in Ref. [32]). They observed also 
an appearance of a void mode as a penetration of a finger- 
shaped (in a vertical section) dusty free region through 
the side boundary of the gas discharge (see Fig. 5c in 
Ref. [32]). A sharp end of the finger and its fast travel 
across the gas discharge volume count in favor of nega- 
tive surface tension of the interface at its end. It means 
that an account for non-spherical form of the void is nec- 
essary for any theoretical model of the void. 
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